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ABSTRACT 



It has been shown that a flux-weighted full width at half maximum timescale of a microlensing event 
can be used in an unbiased estimator of the optical depth. For the first time, this allows a physical 
parameter to be easily estimated from pixel microlensing data. We derive analytic expressions for the 
observed rate of pixel lensing events as a function of the full width at half maximum timescale. This 
contrasts work in the literature which express rates in terms of an "event duration" or Einstein time, 
which require knowledge of the magnification, which is difficult to determine in a pixel event. The full 
width at half maximum is the most directly measured timescale. We apply these results to possible pixel 
lensing surveys, using HST for M87, and CFHT for M31. We predict M87 microlensing rates for the 
HST Advanced Camera and for NGST, and demonstrate that one will be able to probe the stellar IMF. 
Next, we describe a new method by which a crude measurement of the magnification can be made in the 
regime of magnifications A ~ 10 — 100. This in turn gives a crude measurement of the Einstein time. 
This program requires good photometry and sampling in the low magnification tails of an event, but is 
feasible with today's technology. 

Subject headings: cosmology: observations — dark matter — galaxies :halos — gravitational lensing 



The nature of the dark matter in the universe is unknown. The only direct detections to date have been of massive 
baryonic halo objects, MACHOs (Paczyhski 1986), via gravitational microlensing of LMC and SMC stars. The impli- 
cations of the microlensing experiments are controversial, since of the ~ 20 events detected, the 3 events with distance 
determinations are attributable to star-star lensing (2, and the only SMC events, in the SMC, 1 in the LMC). Since 
the SMC is very extended along the line of sight, one must wait for more data to draw any definitive conclusions about 
MACHO location. In the case of the LMC, the distribution of timescales for the events suggests that we should be seeing 
some halo events. If this is the case, up to half of the halo dark matter consists of MACHOs of mean mass ~ O.5M . 

To improve the situation, one must measure many more microlensing events. This can be accomplished by constructing 
much larger cameras and using large telescopes. However there is a complementary strategy that allows us to optimize 
the microlensing yields and that can potentially increase by a factor of two or more the number of detectable MACHO 
events. This is the technique of pixel microlensing (Crotts 1992; Baillon et al. 1993). 

Pixel microlensing works as follows. Individual stars are not resolved. Nevertheless a microlensing event is detectable 
as a pixel brightening. One can observe with many stars per pixel and thereby allow pixel microlensing to provide a 
powerful means of accumulating events. In fact, pixel microlensing allows us to use gravitational microlensing to learn 
about galactic halos even when individual stars are not resolved. The high magnifications required have made it difficult 
to determine the physical parameters of the source-lens system. For example, by observing high magnification events it 
has not hitherto been possible to extract such parameters of the source-lens system as the optical depth. Recently, it has 
been shown (Gondolo 1999) that using the flux-weighted full width at half maximum of the pixel microlensing lightcurve, 
one can extract the optical depth. This motivates us to further explore the theory of lensing in terms of this and related 
timescales. 

We first define some standard notation. The source star being lensed lies a distance L from the observer, and the lens 
itself lies at a distance xL, thus < x < 1. The Einstein radius is the radius of the Einstein ring of a massive object of 
mass M. The Einstein radius and its projected angular size are given by 



Considering the Einstein radius as a function of x, a tube from the source to the observer is traced out. It is widest 
at x = 1/2, and its width goes to zero as yfx or \/l — x as x goes to zero or one, respectively. This is the so-called 
"microlensing tube". A lens that passes inside its own microlensing tube produces a lensing event. 
The Einstein time is the time that a lens would take to cross the full Einstein disk if it were to pass through the center, 
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where v is the lens velocity transverse to the line of sight. The flux increase in a microlensing event is given by the well 
known formula (Einstein 1936) 

=*(„) = /(„»(*)), f(x)= 2 + * -1. (3) 
f V x(4 + x) 

where u(t) is the perpendicular distance from the the lens to the line of sight, in units of the Einstein radius R E . We 
assume linear motion for the lens, and a minimum impact parameter 0, giving u 2 {t) — 2 + (2At/t E ) 2 , with At = t — 1 . 
It is useful to note that when S(u) w l/u. This is the so-called high magnification regime. 

The optical depth to microlensing is a commonly used quantity. It is defined as the number of lenses that are within 
the microlensing tube for a given source. With this definition, the optical depth is clearly given by 



r = I ^R% {x )Ldx = ^^ I p(x)x(l-x)dx. (4) 



Measurement of the optical depth in a microlensing survey provides a weighted column density of lenses, and can be used 
to estimate the total mass in lenses in galactic halos. 

It is straightforward to show (Gondolo 1999) that the full-width at half maximum of a microlensing event is given by 
ifwhm = t E w(0), with 

w(/3) = ^2f(f(0 2 ))-0 2 . (5) 



We find the useful limiting values w(0 <C 1) = 0\/3 and w(0 3> 1) = 0\J \[2 — 1. For classical lensing, where the threshold 
is usually taken to be the Einstein radius, we will need w(l) — 0.9194. The optical depth estimator involves determining 
the flux-weighted timescale t F — Ai^ max ff w h m = Ft E 5(0)w(0). In the regime of high magnifications, required by pixel 
microlensing experiments, -C 1, giving t F w Ft E V3. We use flux to indicate a (dimensionless) number of photons per 
observation. This is simply the number of photons collected in one image of the target. The excess flux is then the number 
of photons in one image above the background. 

2. THE MICROLENSING RATE 

We now would like to derive an expression for the differential microlensing rate in terms of the timescales that can be 
easily measured from pixel lensing data. We have two choices, ifwhm and t F . Recall that it is very difficult to determine 
the Einstein time for a pixel event. For a rate determination alone, the timescale tf w h m is usually more useful than the 
timescale t F . It is precisely the full- width at half maximum timescale of an event that determines if it can be well sampled 
by a given observing strategy. The overall flux comes in to the calculation in determining if a given event can be detected 
above some noise level. If we were to compute the event rates as a function of t F , the sample strategy and photometry 
are tied together in a more complicated way. We thus advocate quoting the event rate distribution in tf w hm rather than 
t F , even though it is the timescale t F which is useful in determining the optical depth from a pixel lensing dataset. The 
event rate is then denoted dT/dtfwhm- This is the differential rate of microlensing events at a given timescale if w hm- 

We modify the rate derivation of Griest (1991) to accommodate the new timescales. The differential rate at which 
lenses enter the microlensing tube is given by 

dT = n(x)f(v- v t )v r d 3 v(2R E d0)(Ldx) = f{v)v 2 r R E Ldv x dv r d/3 dadx. (6) 

Here Vt is the velocity of the microlensing tube in the rest frame of the galactic halos. The tube can move because the 
source and observer arc in general moving in the galactic rest frame. We define v r as the transverse velocity of the lens 
and a as the angle the transverse velocity makes with the velocity of the microlensing tube vt- We quote the impact 
parameter in units of the Einstein radius. We will integrate from zero to some maximum value, which means we 
need the factor of two. We assume that the distribution function f(v) is Maxwcllian, with dispersion velocity v c . The 
distribution function in the microlensing tube frame is 

/(«)=(^)5^«P[-^-W«?]- (7) 

We see that (v — v t ) 2 — v 2 + v 2 + v 2 — 2v r vt cos a, and we can do the integrals in v x and a easily, giving 

dT = Ap( - x ^ LVc Re ( x ) v 2 e - v2 {x) h (2r)v)dv d0 dx, (8) 

with v — v r /v c and r\ = v t /v c . Here, Iq(x) is a modified Bessel function of the first kind. We further define Iq{x) = 
e~ x I (x), thus Iq(x » 1) w l/\/27rx. Previously, this expression has been used to get the rate distribution with Einstein 
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time (Alcock et al. 1995). The variable (3 is trivially integrated to some threshold (3t, and we change variables from v to 
t E . This gives the following, 
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To compute the rate distribution in if w hm, we instead change variables from (3 to if w hm, 2R E d(3 = vv c dt[ w \ lm / 'w' (/?) . The 



transformation is inverted to obtain (3 
be determined later, is then 



w 1 [vv c t[ w h m / (2R E )]. The final theoretical rate, subject to some threshold to 
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If we instead had changed variables from v to if w hm, we would have obtained a double integral in [3 and x, with identical 
functional form in the limit of small (3. We choose to integrate in v and x because it is perhaps more suggestive to 
integrate over the velocity distribution of lenses, then along the line of sight. 

We now must account for the fact that not all events are detectable. The first difficulty is that there must be some 
minimum of additional flux detected so that the event is distinguishable from noise. The simplest way to do this is the 
so-called "peak threshold" . This is the prescription that an event must have some fixed maximum amplification for it 
to be detectable. Effectively, we demand that the signal to noise of the measurement at peak magnification be higher 
than some threshold. We assume that electrons in the CCD are collected at a rate 5ol0~°' 4m , where m is the apparent 
magnitude in some band, and we assume an integration time t- ln t. The signal to noise of the highest magnification sample 
is thus 

Q = <5ma x y^^l0- O - 4(m ^ /2) , (11) 

where m is the apparent magnitude of the unresolved lensed star and p, is the apparent magnitude of the unlensed galaxy 
light on the pixel. We see that a given signal to noise requires some minimum value of 5, which corresponds to a maximum 
impact parameter, (3t- This gives an upper bound on v, as we will see. Fixing all variables except v and (3, we find that 
v is monotonic in f3, thus we arrive at an upper bound on v. 



^c^fwhm 

2R E {x) 



v = w{(3) 



2R E (x) 

111 



w (13 T ) ■ 



(12) 



To improve the peak threshold trigger, we can use the signal to noise in the three highest points of the lightcurve. We 
will assume that the middle sample is on the peak, and that the samples are a time Afti nt apart. The significance of the 
value TV is that the observing program is engaged in on-source integration a fraction 1/7V of the time. The side samples 
are then taken at times Aft- m t off of the peak. The signal to noise as measured this way is given by 



Q 2 = 5 ii„tlO- - 8(m -^ 2) 
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This expression for Q 2 is monotonic in (3, and can be inverted to determine a maximum impact parameter as before. This 
maximum impact parameter then gives the maximum velocity. 

A more sophisticated way to account for the detection of an event is the so-called "event threshold" (Gould 1995). The 
signal to noise in one sample is given by 

= y^t^io- - 4 ^/ 2 ). (w) 

The total signal to noise for the event is thus given by 



(15) 



In converting the sum to an integral, we need to include the factor Af for the observation spacing, since the telescope does 
not integrate continuously. This factor is neglected in previous work (Gould 1995), resulting in overestimates of rates by 
a factor of three or more. Following Gould, we find 
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The correction function ( can be expressed in terms of complete elliptic integrals with parameter m = 4/(4 + (3 2 ). 
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where K and E are complete elliptic integrals of the first and second kinds, respectively. Note that some authors define 
K and E as functions of the modulus k, with to = k 2 . The limiting behavior of this function is given by £(/3 <C 1) = 1 and 
((/3 > 1) = (5/4)/T 6 . For (3 < 0.2, the leading order expansion is accurate to 1%: C(/3) = l-2/3/7r[7r/4-4 + ln(8//3)]. We 
find that when fixing all variables except v, Q 2 {v) is monotonically decreasing. Thus if we impose a threshold signal to 
noise for detection, Qo, we find that w^ ax = Q~ 1 (Qa)- As a final point, in the event threshold case Q 2 is also monotonic 
in P, all other variables being fixed. This means that we can translate the maximum velocity into a maximum impact 
parameter, 

PT = W OD • ( 18 ) 



2R E {x) 

In practice it is quite difficult to implement the event threshold trigger because of non-poissonian errors in the background 
estimation (Gondolo 1998). Using this formalism is more appropriate for space-based observations, which have much more 
easily characterized errors. 

We must now consider the effects of the finite sizes of the source stars. The microlensing formulas assume that the lens 
and the source are pointlike, but this is obviously not the case. The source star has a non-zero radius. The nature of the 
lens is a priori unknown, but the finite lens radius is unimportant for most possibilities. The maximum amplification of 
an extended source assumed to be uniform can be determined by integrating the magnification over the disk of the star. 
The result is 

A^-w-tflgJin, (i 9) 

where the 0's are projected angular sizes of the Einstein ring and the source star, respectively. Following Gould (1995), 
we take the radius of a star of Mj =1.0, 0.5, 0.0,... -3.0 to be R S /R Q = 7, 9, 13, 17, 23, 33, 46, 64, and 90. Such bright red 
stars are predominantly giants. The maximum amplification of a given star is a monotonically decreasing function of x. 
In the peak threshold case, we require that an event has a minimum amplification A m ; n , which corresponds to maximum 
value of x where Anin = ^4max , 
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This is easily accounted for the the rate integral. 

In the event threshold case, we could take the approach that there is a minimum amplification as in the peak threshold 
case. However, we choose to take a very conservative approach. We demand that for an event to be detectable, the lens 
must not pass in front of the disk of the star. This imposes a minimum value of the impact parameter /3 m ; n = s /9e- 
Fixing all other variables, this imposes a minimum velocity, 

-™^IJ0A. (21) 
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In summary, the peak threshold formalism imposes a maximum value of x and a maximum value of v on Eq. [l^, while 
the event threshold imposes both a minimum and maximum value of v on Eq. [l^. 

To get a better feel for the meaning of these results, we compare the bounds on the velocity integral for the peak 
threshold and event threshold triggers in the regime where we require high magnifications. In other words, we can assume 
that j3 <C 1 and 5 = 1/(3. We will assume that we want Q > Qq. Thus, in the peak threshold case, we find 

fr. = 10 -0.4(m- M /2) and ^ = 2 ^ Re(x) V^t ^ 10 -0.4( m - M /2) (22) 

Similarly, in the event threshold case we find 
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For comparison, we form the ratio of these maximum velocities 
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In the limit where the event is well sampled, if w hm 3> A/tint i h is clear that the event threshold prescription allows the 
detection of many more events because the experiment is sensitive to higher velocity lenses. 

Now that we have derived the microlensing rate for a single star with magnitude m which falls on a pixel of magnitude 
H, we need to integrate over the luminosity function of the pixel in order to get the total rate of pixel lensing events. We 
assume that we know the normalized luminosity function of the pixel, 4>(m), in units of stars per unit magnitude. The 
total microlensing rate is then simply 

dT f°° dT 
— =/ dm4>[m)— (rn). (25) 
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For clarity, we now write the total microlensing rate in both the peak threshold and event threshold formalisms, 



dTP 
difwh 



n 



2p Lv 2 
M 



2p Q Lv ( 
M 



dm <fi(m) 



t (m) 



dxg(x) x 



f 

Jo 

/oo />1 
dmcf>{m) I 
-co JO 



Am) 



dvv 3 e 



b>-vWri [2ri(x)v] 



(^c^fwhm 
2R E (x) 



dxg{x) x 

nax( m ) 

v m i n (m) 



(26) 



(27) 



dvv^-^-^Io^ix^} 



( v c U 



whm 



\2R E (a 



These expressions are triple integrals. With a reasonable choice of grid, logarithmic in x or 1 — x, linear in v and m, they 
can be evaluated to 1% accuracy without much computational difficulty. 

For completeness, we would also like an expression for the differential microlensing rate in terms of the flux-weighted 
timescale tp. We return to the differential relation for the rate, and we include the dependence on the luminosity function 
explicitly, 

dT = Ap( - x ^ LVc R E ( x )cp{m)v 2 e-^-^ 2 I [2r){x)v]dv d/3 dx dm. (28) 

We will make a change of variables from m to tp. We find that tp = So*intl 0~ OAm 2R E (x)/ (v c v)w(0)5(0), and furthermore 
that dm, = — 2.5 dtp/(tp In f 0). With the peak threshold triggers, where there is a clear threshold impact parameter 0p, 
we can easily find an expression for the rate, 
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When a high magnification is required, 0t 1 and w(0)S(0) w V3, giving a formula involving one fewer integral, 
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(30) 



As stated before, the timescale t E can be used in an unbiased estimator of the optical depth, as shown by Gondolo 
(1999). However, it is more difficult to interpret the experimental meaning of tp, and it is probably more useful to use 
the timescale if w hm in discussing pixel microlensing rates. Of course, the physical meaning of tp is more clear: in the 
high magnification regime, it is just proportional to the absolute flux of the source star times the Einstein time, with no 
dependence on the impact parameter. 

3. MEASURING THE EINSTEIN TIME 

The flux-weighted full- width at half maximum timescale of an event can be used to determine the optical depth. We are 
still interested, however, in the physical parameters of the source-lens system, encoded in the Einstein time. The difficulty 
in pixel microlensing is that the shape of the lightcurve does not easily allow the determination of the Einstein time at 
high magnifications. In fact, shape alone is not a good indicator of the Einstein time even with small magnifications 
(Gould 1996; Wozniak & Paczyhski 1997). In microlensing with resolved stars, there is extra information in that the flux 
of the source star is measurable. 

With the estimator of Gondolo in mind, we have a situation different than that of previous work. We can measure the 
optical depth to microlensing with no knowledge of the Einstein time. This already gives us a great deal of information 
about the MACHO population in a galaxy. With this information in hand, we can attempt to determine the Einstein 
times for some subset of the data which has high signal to noise. 

In the high magnification regime, we see that 5 w 1/u, giving the shape of the lightcurve 
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where B is the background and F is the flux of the unlensed star. At high magnifications, we see that even several times 
ifwhm off of the peak, the magnification is well characterized by 5 w 1/u. In this regime we can only measure the quantities 
F//3 and f3t E « if w hm/\/3. 

In order to measure the three parameters F, 0, and t F separately, we need to observe the deviations from 6 — 1/u. 
Gould (1996) has shown that this is possible with pixel lensing data with signal to noise Q > 50. To next order in u, the 
lightcurve is 
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This form contains both F/ and F0, breaking the degeneracy. Unfortunately, the F(3 term is quite small, requiring a 
high signal to noise to obtain a robust fit, as Gould describes. 

We describe a new method for measuring the three degenerate parameters in pixel lensing data using another timescale 
that can in principle be directly measured from pixel microlensing data. This is the standard deviation (in time) of 6[u(t)]. 
We define this statistical timescale as follows, 
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We can evaluate £ in terms of complete elliptic integrals, 
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where K(m) and E(m) are the complete elliptic integrals of the first and second kinds as functions of the parameter 

m = 4/(4 + 2 ). We find the limiting behavior £(/3 > 1) = 0/2, £(/3 «1) = [3 ln(8//?) - 6]~ 1/2 , and £(1)=0.778 for 
events tangent to the Einstein ring. 

We now see the utility of t a . At high magnifications, corresponding to small 0, the functional forms of if w hm and t a 
are quite different. If we form the ratio icr/ifwhm, the Einstein times cancel, and we are left purely with a function of 
0. Thus, if we could measure both if w hm and t a for an event, we could measure 0, and thus the magnification and the 
Einstein time. 

We find that t a S> ifwhm in the high magnification regime. In other words, the tails of the lightcurve are very wide. 
This raises the question of how well t a can be measured. It turns out that t a cannot be measured without extending to 
hundreds of times if w hm in the monitoring program. However, we can modify this timescale to make it useful. For a given 
event, we assume that the lightcurve is measured in the interval (— riif w h m , Ji'ifwhm)- We choose to truncate the integrals 
in Eq. ([33) at these multiples of if w hm, for each event. In other words, we measure if w h m for an event, then we compute 
the second moment, truncating the integrals. This defines a new timescale, which we denote the truncated statistical 
timescale, 
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with x n = 2nw(0), f n = [(4 + 2 + x 2 n )/{0 2 + x 2 n )Y 



and 



tan (x n /0). Note that F and E are elliptic integrals 



of parameter to = 4/(4 + 2 ) as before, and that approaches £ for large n. For the full event, with data in the range 
(— nffwhmj ^'ifwhm), we can use the previous formula, and simply duplicate each term in the numerator and denominator 
replacing n with n' . 

The real focus of this method is the following. In pixel microlensing, it is easy to measure the width of the peak, if w hm- 
Events with the same tf w hm have drastically different magnifications. Our method attempts to measure the width of the 
tails of the event, t a , which breaks the degeneracy seen in if w hm- In effect, if w hm sets the overall timescale of the event, 
while t„ splits the degeneracy in events of a given timescale. It is the comparison of the two widths, the scaling of one in 
terms of the other, which allows us to determine the physical parameters of the source-lens system. 

The truncated statistical timescale offers another piece of information about a microlensing event. Because it is only a 
moment of the lightcurve, we can hope to be able to measure it with some accuracy. To illustrate the application of this 
method, we form the appropriate ratio, t^n /ifwhm, and plot versus for n = 10, 20, oo in figure [l| Note that n = oo is 
the idealized case, giving £ CT Afwhm- For n — 10, the ratio varies by a factor of more than three between magnifications of 
a few and magnifications of a hundred. A reasonably good measurement of the ratio can thus give a crude measurement 
of the true magnification. 

4. OBSERVATIONAL POSSIBILITIES 

To illustrate these results we specialize to the case of M87, a giant elliptical in the Virgo cluster. We take the distance 
modulus to Virgo to be D = 31. We use the /-band luminosity function of Terndrup et al., measured for the bulge of the 
Milky Way. We interpolate linearly in (Mi, log N) as follows: (-3.5, 0.0), (-0.6, 1.0), (0.0, 1.5), (0.5, 1.3), (1.25, 1.5), (1.5, 
1.7), (2.5, 2.1), (3.12, 2.65), (4.41, 2.95), (7.35, 3.32), (10.0, 3.9). At the distance modulus of 31, these correspond to to/ 
in the range mi = 27.5 — 41. 

Pixel lensing observations of M87 require HST or the future NGST. For WFPC2 on HST, we find that S = 4.5 x 10 9 
s _1 . As proposed by Gould (1995), we consider using WFPC2, centering the PC chip on the center of the galaxy. With a 
program of dithering, we can achieve a PSF of 1 WFPC2 pixel, which is 0.1" on a side. For simplicity we assume that we 
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can obtain similar resolution on the PC chip, taking the PSF to be 4 PC pixels, each 0.0455" on a side. Since each of the 
three WFC and the single PC chips is 800 pixels on a side, we effectively have 2.08 x 10 6 pixels. M87 is at a relatively low 
declination, which means that in an HST orbit of 96 minutes, we can obtain about tint =1800 s of on-target integration. 
Notice that this implies that Af > 3 for HST observations of M87, since the telescope can only integrate for one third 
of each orbit. This fact is neglected in previous work (Gould 1995). We also compute rates for the Advanced Camera 
(AC) to be installed on HST in the near future. This CCD has 4096 0.05" pixels on a side, and has an efficiency in the 
/-band approximately seven times as large as WFPC2. We take Sq = 3.15 x 10 10 s _1 for the Advanced camera. Finally, 
we compute rates for the proposed Next Generation Space Telescope. We assume a CCD that is 8192 0.03" pixels on a 
side with comparable efficiency to the Advanced camera, and an aperture of 7.2 meters, three times as large as the HST. 
This gives So = 2.84 x 10 11 s" 1 for the NGST. 

We assume that the Milky way halo is an isothermal sphere with a core of 5 kpc and a dispersion velocity of 220 km 
s _1 . The density is normalized to be 0.3 GeV cm -3 at the solar circle. We truncate the halo at a radius of 200 kpc. We 
assume that M87 has a very similar halo, though ten times as massive and with a dispersion velocity of 695 km s . We 
take the core and truncation radii to be 5 kpc and 200 kpc respectively. 

We now compare the rates of observable events in several cases. First, we investigate the peak threshold trigger. In this 
case the rate at which observations are taken only determines the range of timescales of observable microlensing events. 
The sample rate determines the smallest timescale, and the total duration of monitoring determines the largest timescale. 
The rate distributions for M = 0.01, 0.1, 1.0 are plotted for WFPC2, AC, and NGST in figure |. We have assumed that 
Qo = 7 is required for detection. Notice that most observable events are quite short. We also plot the event rates for the 
Advanced Camera with Qo = 42, appropriate for attempting to determine the Einstein time of individual events. 

In figure [| we plot event rates using the improved peak threshold trigger, assuming that we make one observation per day. 
For comparison, we present the rates for the Advanced Camera assuming that we make four equally spaced observations 
in one day. Using the three highest samples in the lightcurve to define an event, we find significantly improved detection 
rates for timescales longer than the spacing of the observations. 

We plot the same four sets of rates using the event threshold trigger in figure ^. Again, we notice a significantly 
enhanced rate of observable events with timescales longer than the spacing between observations. 

Microlensing can also occur when the lens itself is also a star. The rates for star-star lensing are computed in much 
the same way. We assume that a typical star has a mass of 0.4M Q (Gould, Bahcall, & Flynn 1996), and that the stars 
are in an isothermal distribution, with a core radius of 2 kpc. The total mass in stars for M87 is taken to be the same 
as the total dark mass of the Milky Way. It is possible that there are a great number of brown dwarfs in stellar systems. 
We also consider the case where the typical star mass is 0.05 M Q , indicating that most of the mass in stars is in brown 
dwarfs. In figure |5|, we plot these event rates in four cases, all with the improved peak threshold trigger. 

The Andromeda Galaxy, M31, has been the subject of several ground-based pixel microlensing surveys (Crotts 1992; 
Crotts & Tomaney 1996; Tomaney & Crotts 1996; Baillon et al. 1993; Ansari et al. 1997). This galaxy is a somewhat 
peculiar target for a microlensing survey. The disk of the galaxy is sharply inclined, only 12.5° from edge-on. Because 
most microlensed stars lie in the disk, there is a definite bias in the distance to the lensed star. There will be a noticeably 
higher optical depth to stars on the far side of the disk. We take the M31 halo to be identical to the Milky Way halo 
previously described, except that it is twice as massive. In figure O, we show contours of constant optical depth to the 
M31 disk. The isophotes are taken from Hodge $ Kennicutt (1982). The change in optical depth from the near to far 
edge is clearly seen. We now compute event rates for this galaxy. We take the distance modulus of D — 24.3. This 
smaller distance modulus allows us to observe pixel microlensing from a ground based observatory. The CFHT 3.6 meter 
telescope on Mauna Kea is ideal for this purpose. It has a new wide-field CCD, the CFHT12k, which has 12288x8192 0.2" 
pixels, for a field of view of 41' x 27'. We take a PSF of 5 x 5 "super-pixels", making a square arc second, conservative 
for CFHT. We assume that Sq = 5.66 x 10 10 , and furthermore that the noise level is twice the photon counting noise, 
a level obtained by the AGAPE collaboration (Ansari et al. 1997). We use a l/-band luminosity function appropriate 
for the disk of a spiral galaxy (Bahcall & Soneira 1980). Furthermore, we truncate the luminosity function at the bright 
end when it implies a mean of one star per super-pixel of the brightness of the pixel. This is a crude way to remove the 
events involving resolved stars. As before, we compute rates for the peak and improved peak threshold triggers when the 
sampling rate is once per day. We will assume 3600 second exposures. We center the CCD on M31, and align its long 
axis along the major axis of the galaxy. In most cases we take Qo = 7, though we also compute the rates for Qo = 42 
with the improved peak threshold trigger. In addition, there are a large number of resolved stars in this program. We 
compute the rates of resolved microlensing events, assuming that a flux increase of 1% is needed. The rates for CFHT 
observations of M31 are shown in figure [?]. 

In figures ||-^| and [7| we plot dr/c?lnif w h m - This is the differential event rate for events of a specific tf w hm- To get the 
total event rate, we must integrate this over the range of ln£f wnm that the sampling strategy is sensitive to. Roughly 
speaking, we can observe events with if w hm between about the twice the sample spacing and the total monitoring duration. 

5. FEASIBILITY OF EINSTEIN TIME MEASUREMENT 

We now discuss the possibility of using the statistical timescale method to determine the Einstein time. We know that 
we will need a high signal to noise for any method to work. We specialize to the case of a pixel lensing experiment that 
may resolve a great number of stars, perhaps even enough for a classical microlensing survey. This is the case, for example, 
when using the new 12288x8192 pixel CCD on the CFHT to observe M31. The pixel scale is 0."2, giving a substantial 
field of view, with a large number of resolved stars in a modest integration time. With this type of data, the rate of high 
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signal to noise events on unresolved stars is quite high, thus we might hope to be able to recover the Einstein times of a 
subset of the pixel events. 

Measuring even the truncated statistical timescale t an is quite difficult. We now investigate the real- world possibility 
of learning anything from this measure. We perform Monte-Carlo simulations of the observed t an versus the true t an to 
determine the errors involved in this process. 

We proceed as follows. For a given lightcurve in a dataset, we can easily measure ifwhmj which also tells us to what 
value ntf w hm the lightcurve is well measured. We can then fit a constant background to the lightcurve and determine t an . 
We have the theoretical curve for t an , which lets us read off the minimum impact parameter 0, which gives us the full 
information on the microlensing event. 

We would like to know what the error in the determination of is. From the lightcurve, we can determine the signal 
to noise in the peak. We can use a Monte-Carlo simulation to generate events at that fixed signal to noise, for different 
values of beta. For each generated lightcurve, we compute t an , and thus we generate a histogram of measured values, 
which allows us to determine the error in the measured t an , thus an error in the measured [3. 

As a concrete case, we take microlensing lightcurves with peak signal to noise Qq = 42, six times the basic detection 
threshold. We fix tf w iim = 3 days, take exposures of 3600 seconds, and vary combinations of the minimum impact 
parameter f3, and the apparent magnitude my that give our chosen peak signal to noise. Observation are taken once per 
day for thirty days on either side of the maximum. In figure ^, we show generated lightcurves for these types of events. 
For 5 = 20 and S = 100, we give example lightcurves with if w h m = 3 days, for peak signal to noise of Qq = 7 and Qq = 42. 
The values of j3 that we choose correspond to 8 = (10, 20, 50, 100) in the peak. The corresponding absolute V-magnitudes 
are My — (1.78,2.53,3.42,4.28). We now perform a Monte-Carlo simulation in which we generate events of these fixed 
parameters and measure directly the ratio £io = io-ioAfwhm- With these parameters, we find that the four values of /3 are 
distinguishable at the la level. In figure ^, we plot the event histograms for these cases. 

For this technique to work, it is important that a good determination of the background is made. We will assume that 
the background is determined to ler by fitting a degenerate lightcurve. The error will be approximately the square root of 
the number of counts, times the factor of two to account for non-poissonian errors, and divided by the square root of the 
number of samples. We find that this error in the determination of the background induces an error in the determination 
of £io of ±0.029, considerably smaller than the other errors. Thus, the error in the background determination is negligible. 

6. DISCUSSION AND CONCLUSIONS 

The flux-weighted full-width at half maximum timescale can be used in pixel microlensing experiments to measure the 
optical depth. This determination does not require the almost degenerate fit to the lightcurve, as was necessary with 
previous methods. The optical depth is a strong indicator of the amount of mass in lenses along a line of sight. With this 
knowledge, we can further investigate the properties of the lenses. 

It is clear from figures that the shape of the timescale distribution of events depends strongly on the mass of the 

lenses. If a suitable population of events can be detected, it will be possible to make a rough estimate of the mass of 
the intervening lenses. This will serve to provide hints as to the nature of the lenses, which is unknown. For example, 
a large population of solar-mass lenses, inferred to be neutron stars, may be difficult to reconcile with the metallicity of 
the galactic gas. Such a population might be assumed to be primordial black holes. A sufficiently large optical depth of 
any mass lens may also be difficult to reconcile with the total baryon density of the universe. In contrast, microlensing is 
a probe of the low-mass end of the stellar luminosity function, and may give some indication of the number and mass of 
brown dwarfs in galactic systems. 

With a measurement of the optical depth in hand, measuring the Einstein times of a subset of the lightcurves provides 
a great deal more information. This is a more discriminatory determination of the masses of the lensing population. Even 
if only a small subset of events can be so analyzed, the added information is very instructive. 

In conclusion, we have expanded the theory of pixel microlensing in light of the result of Gondolo (1999) that a flux- 
weighted timescale can be used to unambiguously determine the optical depth even with moderate quality data. We have 
described how to compute event rates in terms of the new timescales, and applied the techniques to several proposed 
microlensing surveys. Lastly, we have described a new method by which the Einstein time can be determined from high 
quality data. This method is complementary to Gould's (1996) method. Though complicated to implement, it is perhaps 
conceptually simpler. The physical parameters are determined not through a lightcurve fit, but through a comparison of 
the widths of the peak and the tails. 
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Fig. 1. — The dependence of the peak impact parameter on the ratio t CTn /tfwhm = £,n/ w f° r several values of n. This ratio varies significantly 
for j3 between a few tenths and a hundredth. Thus it may be feasible to make a crude measurement of (3 from measurements of t an and tf w h m - 
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Fig. 2. — Pixel lensing rates for the peak threshold trigger on WFPC2, AC, and NGST. We have taken Qq = 7. In the bottom right we 
plot rates in the Advanced Camera with Qo = 42. 




Fig. 3. — Pixel lensing rates for the improved peak threshold trigger on WFPC2, AC, and NGST. The spacing of observations is one day, 
except in the bottom right plot, where it is six hours. 




Fig. 4. — Pixel lcnsing rates for the event threshold trigger on WFPC2, AC, and NGST. The spacing of observations is one day, except in 
the bottom right plot, where it is six hours. 




Fig. 5. — Star-star lensing rates for the improved peak threshold trigger on WFPC2, AC, and NGST. The spacing of observations is one 
day, except in the bottom right plot, where it is six hours. We have plotted both a typical estimate of the average star mass, 0.4 Mq, along 
with an extreme estimate, implying a great number of brown dwarfs, 0.05 Mq. 
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Fig. 6. — Contours of constant optical depth to the disk of M31. These are given in (dimensionless) units of 10 — 6 . The isophotes (thick 
lines) arc taken from Hodge & Kcnnicutt (1983), and are 19, 20, 21, 22, 23, 24, 25 B magnitudes per square arc second. The rectangular 
boxes show the field of view of the CFHT12k CCD, centered on M31 and aligned along both the major and minor axes. 
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Fig. 7.— Pixel lensing rates for M31 using the CFHT 12kx 8k CCD. The top left shows the peak threshold trigger with Q = 7. The 
top right shows the improved peak threshold trigger, also with Qo = 7. The bottom left shows the improved peak threshold trigger with 
Qo = 42, six times larger than the other cases. The bottom right shows the rate of events where the source star is resolved. The spacing of 
observations is one day. 
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Fig. 8. — Sample lightcurves for M31 events in CFHT. The dashed lines indicate the background level and the prescribed detection levels, 
7<r and 42<r, respectively. The solid lines are the true lightcurves, and are not fits. Top left: Qo = 7, S = 20. Top right: Qo = 7, S = fOO. 
Bottom left: Q = 42, 5 = 20. Bottom right: Qo = 42, <5 = 100. 




Fig. 9. — Histograms of measured t CT lo/tf w hm- 111 a ll cases, tf w h m = 3 days, and the peak signal to noise is Qo = 42. The five dashed 
vertical lines correspond, from left to right, to the true values when <5 = 10, 20, 50, 100, oo. The four points and error bars indicate the mean 
measured values and errors. 



